A New Criterion for Crack Formation in Disordered Materials 



O 

o 
o 

O 

Q 



o 

a 

■4— » 



^3 

o 
o 



> 



(N 

o 
o 

-I— > 
c3 



i 

G 

O 

o 



X 



Peter F. Arndt and Thomas Nattermann 
Institut fur Theoretische Physik, Universitdt zu Koln, Ziilpicher Str. 77, D-50937 Koln, Germany 

(February 1, 2008) 

Crack formation is conventionally described as a nucleation phenomenon despite the fact that the 
temperatures necessary to overcome the nucleation barrier are far too high. In this paper we consider 
the possibility that cracks are created due to the presence of frozen disorder (e.g. heterogeneities or 
frozen dislocations). In particular we calculate the probability for the occurrence of a critical crack 
in a quasi two-dimensional disordered elastic system. It turns out that this probability takes the 
form of an Arrhenius law (as for thermal nucleation) but with the temperature T replaced by an 
effective disorder temperature T e s which depends on the strength of the disorder. The extension of 
these results to d — 3 dimensions is briefly discussed. 



PACS: 62.20.Mk, 61.43.-j 

I. INTRODUCTION 



Cracks are one of the most important sources for the 
failure of solidsa. Despite continuous efforts for more 
than a century aJull understanding of fracture has not 
yet been reachedtm. A simple but very appealing pi6- 
ture for the occurrence of cracks goes back to GriffithM. 
Griffith describes the emergence of cracks as a nucleation 
phenomenon: To open a crack in a thin plate, atomic 
bonds have to be broken and two new surfaces have to 
be created. For a crack of linear size a, this costs an en- 
ergy of the order a. Simultaneously, the potential energy 
of the plate under external load is reduced due to the 
stress relaxation around the crack. This decreases the 
energy by an amount of order a 2 . Thus, the total crack 
energy as a function of a increases for small a linearly and 
reaches a maximum at a = a c before it decreases quadrat- 
ically. Cracks of length a < a c are stable whereas those 
with a > a c are unstable. However, contrary to conven- 
tional nucleation phenomena the typical energy barriers 
for crack propagation in a perfect solid under realistic 
stresses are much too high to be overcome by thermal 
fluctuations. Instead, the pre-existence of micro-cracks 
(or pre-weakened bonds) on scales a < a c is usually tac- 
itly assumed. These will then grow under external load. 
It seems to be reasonable to consider micro-cracks as well 
as other heterogeneities in the material as a kind of frozen 
disorder only amenable to a statistical treatment. 

The propagation of supercritical cracks in an inhomo- 
geneouSj-material has been the subject of a number of 
articlestffl which have attempted to explain the roughness 
of crack fronts found experimentallya. Unfortunately a 
convincing explanation of the experimental data is still 
lacking. The other aspect, the occurrence of a critical 
crack in the first place has been, to the best of our knowl- 
edge, not yet considered. It is the aim of the present 
paper to address this point by calculating the probabil- 
ity of the occurrence of a critical crack in systems which 
includes various types of disorder. In particular, we will 
consider randomness in the atomic bond strength as well 
as randomly distributed impurities (or other kinds of het- 
erogeneities) and frozen dislocations. It should however 



be mentioned, that our considerations are not restricted 
to crystalline material. The main ingredients of our the- 
ory is isotropic elasticity, which also applies to amor- 
phous materials, supplemented by randomly distributed 
disorder. The latter can also include mesoscopic hetero- 
geneities and microcracks which occur during fabrication. 
Although the various sources of disorder conceivable may 
differ considerably in their local properties, the most im- 
portant aspect from the statistical point of view (which 
is adopted in this paper) is the spatial decay of the stress 
fields they create. It is in this sense that the three types 
of disorder considered below are generic. 

The main body of the paper is related to crack cre- 
ation in a thin plate of infinite extension. Some results 
can however be easily extended to bulk materials. The 
detailed investigation of cracks in slabs of finite width as 
well as those in three dimensionaLsystems will be pre- 
sented in forthcoming publications!!! 



II. CRACKS IN AN INFINITE 
TWO-DIMENSIONAL SAMPLE 

In this paper we consider a single planar crack ex- 
tended parallel to the ir-direction in a 2-dimensional elas- 
tically isotropic system of infinite extension. The 2- 
dimensional situation can be realized e.g. by a thin plate 
of width h in the so-called plane stress configuratior£3. 
The Lame coefficients A and /2 of the 2-dimensional sys- 
tem are then related to the Lame coefficients A and fj, of 
the bulk by A = 2A/u/i/(A+2/z), ji = fih. The coordinates 
of the crack are given by 



a < x < a , y = . 



(2.1) 



Such planar crack appears for instance in experiments 
with preweakened bondsEil. 

In two dimensions only mode I and mode II cracks 
occur characterized by the only non-zero components 
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ply the dislocation theory of cracks 
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described by virtual lattice planes filling its interior such 
that there is no free crack surface. The lattice planes 
terminate in crack dislocations with Burgers vector 
The whole crack is then given by a collection of disloca- 
tions (and antidislocations) ba^ at positions r a = (x a , 0) 

Crack dislocations interact with the external stress a. 

y 

and the disorder made up of impurities and frozen dis- 
locations. For the further discussion, it turns out to be 
convenient to introduce a 2-dimensional dislocation den- 
sity b(r) = ^2 a b a S(r — r Q ) . The actual distribution 
of the crack dislocation will be determined later from 
a minimum condition for the free energy for given ex- 
ternal stresses and crystal imperfections. It should be 
mentioned however, that the crack description by dislo- 
cations is not essential for the final results. We could 
also have used more traditional elasticity theory com- 
bined with the appropriate boundary conditions on the 
crack surface. In this sense also amorphous materials 
are included (but there are no frozen dislocations in this 
case). 

The interaction between the external stress a!~ and a 
dislocation with Burgers vector b is given by the Pcach- 
Kohler forceEl fa — — e a b m ai^ , where e« denotes the 
total antisymmetric unit tensor. With the help of this 
relation one obtains for the total contribution of a$ on 
the crack dislocation energy 



E 



-(e) 



-e /'Vi 6 (c) 



fW / dxxbtXx). 



(2.2) 

Thus, in mode I (II) only the y (x) component of b( c )(a;) 
contributes to E^- c \ Since E^ is the only macroscopic 
term which favors the existence of crack dislocations, it is 
clear that in mode I (II) only those with b( c ) (r) parallel 
to the yU$-)axis will occur. This will be used in the 
followingta. 

The stress field er^ generated b 
to the Airy stress function x(r) b 
where x is a solution of 



by-dis 



slocations is related 



(V 2 ) 2 X (r)=Fe^-6 i (r). 



(2.3) 



Here Y — 4/2(A + /2)/(2/2 + A) denotes the Young modulus 
i n tw o dimensions. In an infinite system the solution of 
(2.3) is given byLLii 



X (r)=Y /dV^r-r'K^r'). 



(2.4) 



Here g(r) = r 2 (log |r| + C)/(8ir), C is an arbitrary con- 
stant. The elastic energy of the dislocations is then given 
by 



(2.5) 
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J d 2 r J dV e ii e H 6 i (r)6,(r / )flS i e fcff (r - r') 



(a) 




FIG. 1. Lattice planes terminating in crack dislocations 
in mode I (a) and mode II (b) cracks. The dislocation vectors 
are perpendicular to the lattice planes. The arrows denote 
the direction of the external forces / . Changing /' c - ) crack 
dislocations in mode I and II climb and glide, respectively. 



where we used the relation 
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A 



2/2 4/2(A + /2) 



(2.6) 



The elastic energy of the crack dislocations is then given 
by 



E^ = - — 
8tt 



d.r 



dx' b^(x)b^{x')log 



(2.7) 



&( c - ) (x) denotes the crack dislocation density along the x- 
axis and ao denotes a microscopic cut-off (of the order 
of the lattice spacing and in general different for mode 
I and mode II). Eq. (2.7) is valid both for mode I and 
mode II cracks. Note that E^ is always positive. 

So far we did not consider the core contributions of the 
dislocations. In the present context this is replaced by 
the 2-dimensional energy density 7 (x) = 70 + 71(21) of 
the crack surface 



E (s 



dx^(x) = 4j a + E[ s) (a) 



(2. 



71 (x) reflects the randomness in the strength of the bonds 
broken. For simplicity we assume Gaussian disorder with 
(71(2;)) = and 



(71(^)71(^0) = ^7 a 8 aa (,x - x') , 



(2.9) 



(...) denotes the average over the disorder and S ao (x) 
a delta function of width ag. In general th e correlation 
length do of the disorder appearing in ( [2.9] ) is different 
from the cut-off appearing in (2.7). Similarly further cor- 
relation lengths could be introduced for the distribution 
of impurities and frozen dislocations to be considered be- 
low. To avoid a too clumsy notation we will however use 
everywhere the length ao as a small scale cut-off but keep- 
ing in mind this complication. Depending on the type of 
material under consideration, ao may vary between the 
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size of an atom in crystalline and the size of a grain in 
granular materials, respectively. The precise value of ciq 
will be of course important if comparison with experi- 
ments is made. Then (^E^\a)J = and we find for the 

variance of e[ s '(a) 



E[ s \a)-E?>(a>) 



00, 



= A 



00! 



(2.10) 



where A/ s ) = 8<57 ao. Clearly Eq. (2.E) is only valid for 
\a — a'\ > ao- 

In the following we add the contributions E^ d ' — 
£;( fd ) -\-E0) from randomly frozen dislocations and impu- 
rities to the energy. E^ is given by Jd 2 r u\j where 

a^p denotes the stress generated by the disorder and uf^ 
the strain field generated from cracks, respectively. 



Using Eq. (2.5) E^ can be written in the form 



E (d) 



Y 

47T 



dxb {c) (x)V(x) 



V(x) = V (id) {x) + V ({) {x) 



(2.11) 



(2.12) 



From (2J3), ( |2.4|) and (|2.6| ) one obtains for the potential 
created from dislocations 



V^(x) = 4nJ d 2 r' e tJ (d k d l9 (r - v')) y=Q bf \v') 

(2.13) 

where k — x,y for mode I, II cracks, respectively. The 
frozen dislocations are assumed to have both random po- 
sitions r Q and directions of their Burgers vectors such 
that (b( fd )(r)) = and 

(bf d) (r)bf V)) = b 2 {m) c m S ao (r - r')S tJ (2.14) 

Here c^) and b^j denote the concentration and the 
strength of the dislocation. 

Impurities (or more macroscopic inclusions) also gen- 
erate a long-range elastic displacement field u^O(r). Re- 
peating the calculation of Eshclbytj for d = 2 dimensions 
one finds for the strain tensor of an impurity located at 
the origin 



(2.15) 



Here d denotes the 2-dimensional volume change due to 
the impurity which can be of either sign. The interaction 
energy between the crack dislocations and the impurities 
of density C(;) (r) takes the form 



£« = ^/d 2 r4 c) (r)c (i) (r) 



(2.16) 



Here we used ( |2.15| ) and 

5 «( r ) = X^ a «( r ~ rQ ) _S ( i ) 



(2.17) 



where the summation is over all impurity si tes r a and 
denotes the impurity concentration. With ( |2.3| )-(2.6) we 
find 



U (i) (.x) = Vl J dVc(r') (d k \og 



(2.18) 



again k — x,y for mode I, II cracks, respectively. 

The total energy E = E^ + + + E^ is a 
functional of the crack dislocation density £>' c ) (x) . Differ- 
entiating the saddle point equation SE /5b^ c \x) — with 
respect to s, we find 



47T . 

~Y 



-a^ + dx'b^(x') + V'(x)=0 (2.19) 

J -a X-X' 

Eq. (J2.19D has the solutioiS 



b^ c) {x) = f dx' f{x,x';a){^a {c) +V'{x')) 

J —a * 

= bi c \x) + b^\x) 



(2.20) 



where 



f(x,x';a) 



1 a 



,2 _ ^2x1/2 1 



(2.21) 



The total energy as a function of the crack length follows 
with the help of (|1|) 



dx 



dx' b (c \x)b (c \x')lo£ 



«o 



+2 / dxj(x) 



(2.22) 



where b^ c \x) given by (2.20). 

For vanishi ng di so rder b 00 (x) — > bffl (x, a) for which we 
obtain from d^20| ), ( f2~2ll ) and (|aT| ) 



6 c) (x,a) 



Y (a 2 -x 2 y/ 2 



(2.23) 



which yields in mode I an elliptic crack of maximal height 
2d^a/Y. 

As follows from (2.2C) the total energy ( 2.22j ) can be 
divided into contributions E n which are proportional to 
(<7( e )) 2-n with n = 0,1,2, respectively. Here E\ and 
E2 depend on the disorder. The disorder-independent 
contributions to the energy are given by the Griffith ex- 
pression 



na 2 af ) ( a 

E (a) = 470a =M = 4 7o a 1 - — 

Y \ 2a c 



(2.24) 
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which shows a maximum at a — a c — 2jqY / (ira^) cor- 
responding to an energy barrier £o(a c ) = 2^a c . 

The contributions E\ and Ei depend on the frozen dis- 
order and can be characterized by their mean value and 
varia nce . E\ c an be rewritten using partial integration 
and and (|a|) as 



E 1 



Y 

4-7T 

Y 

47T 



—a J —a 



dxdx' b^^b^ (x')log 
dxV(x)b ( g\x) 



X — X 



Of) 



(2.25) 



For impurities we obtain from ( 2.17 ) and ( 2.18 ) 

R 



V {i) (x)V® (a/)) = n 2 c (i) irlog 



\x — x'\ 



(2.26) 



where R is a cut-off of the order of the syst em size which 
has to be send to infinity. This gives with (Al)— (A.2) for 



the impurity contribution E± to the variance of E\ 



Ef{a)-Ef{o!) ~)=A (i) |a 2 - a ' 2 



(2.27) 



Here A (i) = (7r/2)c (i) (ficr (e) ) 



)^ 2 7o£Z 



For frozen dislocations we get from ( 2.13 ) and ( 2.14 ) 

V^(x)V^(x')) = Q -c (id) b% d) (&f- -{x- x'f) 

(2.28) 



which gives with ( 2.25| ) for the dislocation contribution 



E^ to the variance of E\ 



E[ id \a))-Er>(a') ) = A (fd) (a 2 - a' 2 ) 2 (2.29) 



n(fd) , 



with 

A(fd) = (3/7r)c (fd) (6 ( f d) CT(°)) 2 = 6c (fd) 6 2 fd) 7oF/(7r 2 a c ). It 
is easy to see from ([2.11) and the condition that the is 



crack closed, i.e. J_ dxb^(x) = 0, that there depen- 
dence on R vanishes in ( p. 27 ), (2.29). These equations 
are clearly valid only for \a— a'\ larger than ag. For \a— a'\ 
smaller than the mean distance between the impurities 
or dislocations, respectively, t he st ati stics of E\(a) is not 
longer Gaussian, but Eqs. ( 2.27 ), ( 2.29 ) still give the 
correct order of magnitude of the fluctuations of E\ (a) . 

The mean values as well as the fluctuations of Ei are 
proportional to C(;) and C(f d ) and hence small if the dis- 
order is weak as we will assume in the following. Then 
the average energy of the crack is given by the Griffith 
expression Eo(a). The energy barrier Eo(a c ) = 2joa c is 
typically large and cannot be overcome by thermal fluc- 
tuations. Indeed, for crystalline solids witbla 70 ^ Yclq 
one finds 



E {a c ) < k B T m {Y/af 



T m = 7oao/fcB (2.30) 




FIG. 2. Schematic plot of the Crack energy as a function 
of a in the pure and random system. 

where T m is a characteristic temperature compara- 
ble to, but j-typically bigger than, the solid's melting 
temperaturetj. For relatively large strain Y/a is of the 
order 10 such that the nucleation rate for a supercritical 
crack is of the order ojq exp(— 100T m /T)! uj$ is a micro- 
scopic attempt frequency of the order 10 13 s _1 . In the 
further discussion we will therefore mostly neglect ther- 
mal fluctuations. 

Let us denote the probability that a crack of length a 
has a negative energy by We<o{o). The smallest crack 
one can think of has a length of the order a . Thus a 
crack can only appear if this smallest crack has a neg- 
ative energy, E(ao) < 0. This occurs with the proba- 
bility We<o{clq)- (On this smallest scale even thermal 
fluctuations may help to create a crack as we will see be- 
low). The crack can then only propagate further, if for 
a given disorder configuration the force on the crack tip 
f(a) = —dE/da is positive for all a > ao (we neglect 
effects of inertia). 

Because of its mathematical simplicity we consider 
here first a necessary condition to be fulfilled which is 
E(a) < for all a. Clearly, if E(a) > an (essentially) 
macroscopic energy barrier exist and the crack cannot 
propagate. The probability We<o{o) is given by 



with 



W E <o(a) 



4>{a) 



1 



dxe"* 2 / 2 



47 a[l - a/(2a c )] 



(A (s) a + A (i) a 2 + A (fd) a 4 )i/2 



(2.31) 



(2.32) 



Here we used the fact that E\ 
sian distributed with variance A 



00 a 



(0 



E^ is Gaus- 

4 



A (i) a 2 



a 



Note that this expression for the variance makes sense 
only for a much larger than ao- Among the cracks of 
various lengths a there is one length a c for which the 
probability We<o is minimal. This minimum appears at 
the maximum of 4>(a). One finds a c given by the solution 
of 



4 



0=(1 



3a , . a 2 
) A (s) - -A (i ) 



2a 



2a d A 



(id) 



(2.33) 



fo(a) = 47o(l 



(2.39) 



It is t empt ing to consider We<q{o>c) = We<o (compare 
.31), (2.33)) as the probability for the occurrence of the 
crack. This conclusion will be elucidated further below. 
It is instructive to consider the three different sources of 
disorder separately: 

(i) If only the surface energy is random, then 0( s ) (a) 
vanishes at a = and a = 2a c and has an maxi- 
mum at a c = (2/3)a c . For weak disorder </>( s )(a c ) = 
(f>( s ) is large and since for 0(a) 3> 1, We<o{o) ~ 
-exp(-0 2 /2)/(V2^0) we get 



W, 



1 



£<0 



^(j)^ (a c 



■ exp 



167oa e 
27^7 2 ao 



(2.34) 



(ii) In the case of randomly distributed impurities 0(;) (a) 
has a maximum for vanishing a. For weak disorder and 
a c clq we get 



1 



27r0W(o c ) 



exp 



870 a c 



(2.35) 



(iii) Finally, for frozen dislocations 0(fd)(a) = 4to(1/<z — 
l/(2a c ))/ -y/A(f(j). Hence </>(fd)(a) takes its maximum 
again at a — > where 0(a) diverges. In this case a has 
to be replaced by the minimal lengthscale ao- This gives 



W, 



(I'd) 



1 



_B<0 



27T(/)( fd )(a c ) 



exp 



47r 2 7 a c 



36 2 fd) c (fd )a 2 y 



(2.36) 



It may be worthwhile to mention, that thermal fluc- 
tuations could be at least partially incorporated into the 
present treatment by relaxing the condition E(a) < 
to E(a) < fceTlog^o^) since barriers of order E are 
overcome on time scales of the order t w ojq exp(E/T). 
This replacement changes the numerator of 0(a) in eq. 
Q2.32D from 47 a(l - a/(2a c )) to 47 a(l - a/(2a c )) - 
k^T \og(u)Qt) . Since barriers on scales a < or with 



«t « a -— log(w i) 

4^ m 



(2.37) 



disappear, ao has to be replaced by max(ao,ar) in the 
expressions for We<o (we still assume a c 3> «t). This re- 
placement eff ects e ssentially only the result for the frozen 
dislocations, ( 2.36| ). We emphasize however, that in gen- 
eral the three types of disorder will work in parallel. 

Next we consider the probability to fulfill the sufficient 
condition dE/da <0for0<a<oo,i.e. that the force 
/(a) acting on the crack tip is always positive and the 
crack can propagate forever. We decompose / into a 
deterministic and a stochastic contribution:. 



f(a) 



dE(a) 
da 



/ (a) + / 1 (a) 



(2.38) 



h(a) = f i 1 s) (a) + n i >(a) + fr > (a) 



f(fd), 



(2.40) 



Since /i(a) = —dE/da is Gaussian distributed, the joint 
probability distribution of /i(a) for < a < 00 is also 
Gaussian but in general non-local. Its form can be recon- 
structed from the second moments 

(/i ( %)/i (s V)) =A (s) 5 ao (a-a') 

(/Pto/JV)) =2aA (i) Mo-a') 
(/f)(a)A (fd V)) =4aa'A (fd) 

In the case of random surface tension and randomly dis- 
tributed impurities correlations are local and the joint 
probability distribution of the f's factorizes. In this case 
the probability that the force /(a) on the tip of a crack 
of length a is positive is given by 



W f>0 (a) 



-Ma)/(fl(a)Y 



d.r 



1 



-x 2 /2 



2tt 



(2.41) 



The total probability Wf>o for /(a) > for < a < 00 
is given by the product of all W/>o(a) factors. Here 
we take into account, that the forces are correlated over 
a distance of order ao such that we can decompose the 
crack in pieces of the order ao- It is more convenient to 
consider the logarithm of W/>o 



logW/ >0 = ^logW/ >0 (a) 



(2.42) 



The sum over a = nao is here meant over integer num- 
bers n. Since for a 3> a c the integrals (2.41) are essential 
equal to unity, it is sufficient to restrict the summation 



111 ( pl2; ) to the region < a < a c . Moreover, the sum is 



dominated by the smallest values of Wf>o( a) for which 
we can replace the Gaussian integral (2.41) by the ap- 
proximate expression used above. This gives 



logW) 



>o 



° c da / / 2 (a) 



log 



<A 2 («)> 



/„ 2a \{fl(a)) & 27rfl(a)J 

(2.43) 

In the case of a random surface tension only we obtain 



logWg 



Oc I 

a o \ 307 2 



1 , 6j 
1+ 2 l0g i^ 



(2.44) 



Thus W/>o is essentially of the same form as We<o apart 
from a replacement of the numerical factor in the expo- 
nent (^| is replaced by ^). A similar calculation for the 
case of random impurities gives 



5 



(log — - r) 



a \ c (i )fi 2 y 



«o 



ii og «!T| 

2 327ra 7o y 



(2.45) 



which is again the same result for the exponent as for 
We<o, apart from the logarithmic factor which replaces 
2. 

The third case of randomly distributed dislocations is 
more involved. Here we should indeed take into account 
the long range correlations of the forces f[ id ^. This re- 
quires a more detailed mathematical investigation which 
we leave for further studies. However our experience 
with the two other cases makes it tempting to assume 
that W£% gives essentially the correct expression of the 
probability for the occurrence of a crack. 



III. DISCUSSION AND CONCLUSIONS 

It is interesting to remark that in all cases considered 
above W can be written in the form of an Arrhenius 
law for thermal nucleation with an effective temperature 
determined by the strength of the disorder: 



W « Wo exp - 



/ _ 27 q e \ 



Here 



T 



(s ) _ 6(57 q 



off 



fcB7o 

for random surface tension, 

n 2 c (i) Y 



T 



(i) 



off 



4fc F 



for randomly distributed impurities and 



T 



(fd) ^ 3 a C>k(fd) C (fd) 



off 



27T 2 fc F 



(3.1) 



(3.2) 



(3.3) 



(3.4) 



for frozen dislocations. The present calculation is not 
accurate enough to determine the pre-exponential term 
Wo, which we assume to be of the order one. Relation 
(3.1) can be given a very simple meaning: the probabil- 
ity that a crack of minimal length ao occurs is given by 
cxp(— 27oao/(fcB?cff)). Here 27oao denotes the energy of 
such a crack. The total probability is then the (a c /ao)-th 
power of this elementary probability. 

Below we want to estimate W for two different mate- 
rials. We have to keep in mind, that our calculation was 
strictly two-dimensional. To compare the results with 
real experiments on thin plates we have to consider their 
dependence on the width h of the plate. A necessary con- 
dition for the application of the two-dimensional theory 
is, that the critical crack length a c is much larger than 



h. In the following we will assume, that this condition is 
always fulfilled. 

Since 70 = 70ft-, Y — Yh and aY e ) = o~r e \h are propor- 
tional to h, the nucleation energy 27 a c = 47qY / '(ndf^) 
is also proportional to h. 

Estimating the /i-dependence of T c g we have to make 
sure, that all relevant length scales in the xy-plane (like 
a c ) are much larger than the h. For T e g and which 
were determined by the small scale cut-off ao, this leads 
to the severe restriction h < ao. We note however, that 
this cut-off will be in general larger than the lattice spac- 
ing. It is therefore appropriate to use in these two cases h 
as the small scale cut-off. The behavior on even smaller 
scales is described by three-dimensional physics which we 
will discuss at the end of the section. 

With J7 6^ 2 ha , where £7 denotes the fluctuation 

(s) 

of the three-dimensional surface tension 70, T ff does not 
depend on h. 

Next, C(j) = cn\h is a two-dimensional density propor- 
tional to h (the three-dimensional density c/^ is indepen- 
dent of h). n as a two-dimensional cross section of the 
impurity should be replaced by fl 3 / 2 /h, i.e. we assume 
essentially spherical impurities. Thus 



T 



(i) 



err 



n 3 c (i) Y 
4fc R 



V 2a c 



(3.5) 



is essentially independent of /i(<C a c ). 

For randomly distributed dislocations C( fd ) (as a line 
density) and b^) are unchanged and hence 



CIT 



h %d) c (td)Y 



4ir 2 k B 

In ( |3.5| ) and ( |3.6| ) we replaced a by h. 



(3.6) 



On scales smaller than h the cracks are three- 
dimensional. The results of this paper can however easily 
be extended to penny cracks in d dimensions. The Grif- 
fith energy then takes the form 



E (a) 



(e) 



(3.7) 



where we neglect here and in the following all numerical 
factors which in general depend on the precise crack ge- 
ometry. Eo(a) has a maximum at a c ~ loY/c 2 c y i.e. a c 
is essentially unchanged. 

If the disorder is taken into account, additional con- 
tributions to the energy appear. A randomness in 
the surface tension leads to a fluctuation of the order 
^(E[ s \a)) 2 '^ « 5j 2 (aa ) d ^ 1 . Following the discussion 

below Eq. (2.31), the minimum of We<o{o) follows again 
for a c r-j a c . Rewriting the minimum We<o,<i of We<o{o) 
in the form 



W d « Waq exp 



(3. 
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we obtain for the effective disorder temperature (d) « 

57 2 afV(fcB7o_). , 

Randomly distributed impurities create an additional 
contribution to the stress where 



7 «l 



Y{n d c {i) /a d ) 1 ' 2 



(3.9) 



This expression can be understood as follows: An isp=. 
lated impurity creates in a volume a d an average stresstj 
of the order Yfl d / 2 /a d . With c^a d the number of 
impurities in this volume (the average stress created 
by the impurities is assumed to be already incorpo- 
rated into Eo(a)) the fluctuatio ns o f the stress created 
by the impurities is given by (3.9). Thus we obtain 

((Ef\a)) 2 ^ w n d Yj a d c (i) /a c . The minimal probabil- 
ity follows for d > 2 (contrary to the 2-dimensional case) 
for a c r*> a c . The disorder temperature is now of the or- 
der T®(d) = n d c (i) Y/(4k B ). Clearly, for a c « h, T^(d) 
agrees with the result for the plate, as it should be. 

Similarly, randomly distributed dislocation lines will 
give a fluctuation contribution er( fd ) to the stress of the 
order 



T({d) I 



Yb 



1/2 
(fd) C (fd) 



(3.10) 



Indeed, a single dislocation (line) creates a stress of the 
order Ybo/a. With c^a 2 for the number of dislocation 
lines in the volume a d which are a ssumed to have ran- 
dom orientations one obtains ( 3.10| ). The corresponding 
fluctuation of the energy due to dislocations is therefore 

given (^(E[ id \a)) 2 ' K j « a 2d 7o^» 2 fd )C(fd)/a c - The minimal 
probability follows as for d = 2 from small scales a « ao, 
which results in 



r (fd) 

( eff 



(d) 



,2„<i-2 



y& (fd) c (fd)/^B 



(3.11) 



Clearly, also the case o f mul tiple disorder can be con- 
sidered as it was done in ( 2.31 ). 

Next we have to compare the different probabilities in 
order to decide, which process dominates the formation of 
cracks in thin plates. For systems with a random surface 
tension, and under the condition a c h, the minimal 
probability W arises both in two and in three dimensions 
from cracks of length d c ~ a c . In this case the probability 
for crack formation is given by ( |3.1| ) , ( |3.2| ) . 

For randomly distributed impurities d c (d — 2) ss ao ~ 
h in two dimensions and d c (d = 3) ~ a c in three di- 
mensions. Thus crack formation is dominated here by 
two-dimensional cracks of minimal length ao ~ h. The 



corresponding probability is now given by (3.1), (3.5). 

Finally, for randomly distributed dislocation, crack for- 
mation in both two and three dimensions is controlled by 
the formation of small cracks of size d c ao- Compar- 
ing the corresponding probabilities in two and in three 
dimension (for d — 2 we have to use ao ~ h) we find that 
crack formation is dominat ed b y pen ny cr acks of size ao- 
Its probability is given by (|3.8|) and (3.11) for d = 3. 





Glass 


SiC 


Y 


[10 y Pa] 70 
[Jm- 2 ] 1.0 


400 


7o 


4.0 


random surface energy 


weak disorder: 


J 7 / 7o = 0.1, a = 5 ■ 10- iu m 




rW(d = 2) 


[K] 1087 
[PaV^ 1 ] 5.9-10 30 


4348 
1.36-10 32 


strong disorder: S-y/jo = 0.3, ao = 10 °m 




rW(d = 2) 


[K] 3.9- 10 10 
[Pa 2 m _1 ] 1.65-10 23 


1.57-10 11 
3.77- 10 24 


random impurities 


weak disorder: 


f2 = 2.5-10 _lu m 2 ,c = 8-10 24 m" 




T$(d = 2) 
A ® 


[K] 158.5 
[Pa 2 nr 1 ] 4.07-10 31 


905.8 
6.5-10 32 


strong disorder: £1 = 10 ib m 2 ,c=10 i 'm 3 




A ® 


[K] 1.26-10 8 
[PaV" 1 ] 5.09- 10 28 


7.24- 10 5 
4.15-10 29 


random frozen dislocations 


weak disorder: 


6 (fd) = 5-10~ iu m,c (fd) = 10 i4 m 


~' 2 ,h = 10 _3 m 


Tiff' (Plate) 
ao = 5-10 9 m 


[K] 


1.83-10 19 


A <tQ 


[Pa 2 !^ 1 ] 


1.71-10 30 


strong disorder: &(fd)=5T0 iu m,C(fd) = 10 ib m 


~' z ,h = 10 _3 m 


Tiff' (P^te) 
a = 5-10" 10 m 


[K] 


1.83-10 21 


A (fd) 


[Pa 2 m" 1 ] 


1.71-10 31 



TABLE I. Estimates of the effective temperatures T e g for 
Glass and SiC using (3.2) , (|3.5| ) and ([3.6]). The corresponding 
factors A in (3.12) and ( |3. 13{ ) are also given. The material 
constants are taken from Ref. |^. 



In a macroscopic sample of linear size L regions of dis- 
tance greater than d c can be considered to be essentially 
independent. The total probability for crack formation is 
given by Wd(L/d c ) d where d = 2 for surface-and impurity 
disorder and d = 3 for frozen dislocations .Ej 

Next we consider the probability for crack formation 
in two different materials, one is amorphous (glass) and 
one crystalline (SiC). The corresponding parameters are 
summarized in Table 1. In the case of random surface 
tension and randomly distributed impurities we express 
the probability for crack formation as a function ofh/a 2 ^ 
since the probability can be written in the form 



log(wywb) 



4 7o 2 y h _ A 
7T k B T cS cr ( 2 e) a 



h_ 

2 

(e) 



(3.12) 



For frozen dislocations we have 

log(wywb) : 



o 2 



a & (fd) C (fd) %) 



(3.13) 



'(e) 



The corresponding values for A^ s \ AW and A^ are also 
given in Table [j]. 

A number of comments are in order: 
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(i) The results so far are based on the assumption of 
short range correlations of the disorder (surface energy, 
impurities, frozen dislocations). Experimentally this may 
not be the most important situation. Long range corre- 
lations of the di sord er desc ribed by a power law decay 
on the r.h.s. of (2.9), ( 2.14 ) and ( 2.16 ) would in general 
lead to an increase of the probability for the occurrence 
of a crack. 

(ii) In the case of impurities and frozen dislocations the 
probability W is dominated by the energetics on small 
length scales. Since the stress on the crack tip diverges 
in the linear elasticity theory used throughout the pa- 
per, non-linear effects may be of particular importance 
for small cracks. This could diminish the numerical co- 
efficients in the exponent of W . 

(iii) Defects were assumed everywhere to be frozen. 
In real fatigue experiments often alternating stress is ap- 
plied which leads to an accumulation of dislocations close 
to the crack tip, which makes crack propagation easier. 
This mechanism could also help in a situation, were due 
to disorder fluctuations, the stress is considerably higher 
than in the average (or the surface energy is lower). 

(iv) An interesting question is the relation of our re- 
sults to crack propagation in quasicrystalslla, which take 
an intermediate position between periodic and random 
media. 
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APPENDIX A: RELEVANT INTEGRALS 

= — 7t.t for lad < 1 



h= f'dx'ii-xy/ 2 ^— 

J-l x'-x 



h = t dx'il-x 12 )- 1 / 2 ^- 
J-l x — 



for \x\ < 1 



7T — 7T \X \(X 



(Al) 



\x\(x 2 -iy^ 2 for |rc| > 1 (A2) 
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